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^ ■ Proton momentum distributions of the ^^C, ^^O, "'"Ca, ^^Fe and ^°*Pb nuclei are 
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calculated by a model using the natural orbital representation and the experimental data 

■ 

o 
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' distributions to be obtained using only hole-state natural orbitals (or mean-field single- 

(N ■ 

particle wave functions as a good approximation to them). To demonstrate the model 
two different sets of wave functions were employed and the predictions were compared 



> 

o 

, with the available empirical data and other theoretical results 

o 
in 

0\ • I. INTRODUCTION 



■ The systematic investigations of the nucleon momentum distributions in nuclei extend the scope of the 

P5 . nuclear ground-state theory. Until the mid-seventies more attention in the theory had been paid to the 

, study of quantities such as the binding energy and the nuclear density distribution p{r). This is related 

X, 

H to the ability of the widely used Hartree-Fock theory to describe successfully these quantities, which, 

however, are not very sensitive to the dynamical short-range correlations. The experimental situation 
in recent years concerning the interaction of particles with nuclei at high energies, in particular the 
(p, 2p), (e, e'p) and (e, e') reactions, the nuclear photoeffect, meson absorption by nuclei, inclusive proton 
production in proton- nucleus collisions, and even some phenomena at low energies such as giant multipole 
resonances, makes it possible to study additional quantities. One of them is the nucleon momentum 
distribution n{k) which is specifically related to the processes mentioned above. However, it has 
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been shown |g| that, in principle, it is impossible to describe correctly both momentum and density 
distributions simultaneously in the Hartree-Fock theory. The reason is that the nucleon momentum 
distribution is sensitive to short-range and tensor nucleon-nucleon correlations. It reflects the peculiarities 
of the nucleon-nucleon forces at short distances which are not included in the Hartree-Fock theory. This 
requires a correct simultaneous description of both related distributions p{r) and n{k) in the framework 
of nuclear correlation methods. 

The main characteristic feature of the nucleon momentum distribution obtained by various correlation 
methods [1, 2, 4-23] is the existence of high-momentum components, for momenta k > 2 fm^^, due to 
the presence of short-range and tensor nucleon correlations. This feature of n{k) has been confirmed 
by the experimental data on inclusive and exclusive electron scattering on nuclei (e.g. [1, 2, 24-27]). 
We emphasize also the fact that theoretical results of various correlation methods [p|,p^ 
experimental data for n{k) obtained by the y-scaling analysis [p6| confirm the conclusion ||5|,p8|] that the 
high- momentum behaviour of the nucleon momentum distribution {n{k)/A at fc > 2 fm^^) is similar for 
nuclei with mass number A=2, 3, 4, 12, 16, 40, 56 and for nuclear matter (see p. 139). More precisely, 
the high-momentum tails of n{k) are almost the same for all nuclei with A > 4 and thus ^He is the 
lightest nuclear system that exhibits the correlation effects via the high-momentum components of the 
nucleon momentum distribution. Since the magnitude of the high-momentum tail is proportional to the 
number of particles, this effect is associated with the nuclear interior rather than with the nuclear surface. 
This allows us in the present paper to suggest a practical method to calculate the proton momentum 
distribution for nuclei heavier than "^He (e.g. ^^C, ^^O, '^'^Ca, ^^Fe and ^'~'^Pb) from that one of '^He 
which is already known from the experimental data (or from calculations within correlation methods |^). 
Here we should like to emphasize that though our method has some similarities to that one suggested 
in ||l6| (extended and developed in ||29|-|3l[] ) , in contrast with the previous calculations, the correlation 
effects are extracted from ^i/e rather than from nuclear matter. We should like to mention also that the 
experimental data for n{k) in ^He (which we use in our calculations) as well as for other nuclei (which 
we use for a comparison) are not directly measured but are obtained by means of the y-scaling analysis 
[ p6[ relying on the assumption that the l/q expansion is valid. For this reason we give an additional 
comparison of the data for n{k) in ^i/e j2^ with the theoretical calculations from p9[ ]. 

In general, the knowledge of the momentum distribution for any nucleus is important for calculations 
of cross-sections of various kinds of nuclear reactions. It is known that reliable results for n(k) in sophis- 



as well as 
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ticatcd methods such as the exp{S)-mcthod g, the variational method with state dependent correlations 
[ pO[ , the generator coordinate method (with two generator coordinates |^^) and others are available only 
for light nuclei up to the ^^O nucleus. The local density approximation (LDA) which has been applied 
to derive the spectral function of finite nuclei and to calculate n{k) in ^^O was used as a basis to 
calculate n{k) also in ^°Ca @ and in nuclei with A=16, 40, 48, 90 and 208 Q. 

The model suggested in this work uses the transparency of the single-particle picture being within the 
framework of a given correlation method by means of the natural orbital representation | ]33[ |. The latter 
enables us to specify in a natural way the high-momentum components in the momentum distribution 
which are of the same magnitude for various nuclear systems. The theoretical scheme of the method 
combines the mean-field predictions for the nucleon momentum distribution which are expected to be 
realistic at small k (k < 2 fm^^) with the correlated part of the momentum distribution. In this sense, 
the method in this work has a similarity with that one from |29[ proposed for calculations of the spectral 
function P(k, whose energy integral the nucleon momentum distribution is. We emphasize that in 
our work this is done upon the common ground of the natural orbital representation. The analyses of 
n{k) performed in this work which use essentially the correlations contained in '^He nucleus (in contrast 
with the calculations based on nuclear matter results already mentioned) can diminish the theoretical 
uncertainties on n{k) for medium-heavy nuclei. 



II. THE MODEL 



We start from the natural orbital representation |33| , where the proton momentum distribution nor- 
malized to unity is of the form : 



nlj 



Rnljik) 



(1) 



where Xnij is the natural occupation number for the state with quantum numbers (n, l,j) and 



^(2j + 1)A„,, = Z. 

nlj 



(2) 



The radial part of the natural orbital in the momentum space Rnij{k) is related to the radial part of the 
natural orbital in the coordinate space Rnij (r) by 



R,uj{k) = (2/^)i/2(-*)' / r^ji{kr)Rni,ir)dr, 
Jo 



(3) 
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where ji{kr) is the spherical Bessel function of order I. We call hole-state natural orbitals those natural 
orbitals for which the numbers Xnij are significantly larger than the remaining ones, called particle-state 
natural orbitals [Q. It was shown by the Jastrow correlation method that the high- momentum 
components of the total n{k) caused by short-range correlations are almost completely determined by the 
contributions of the particle-state natural orbitals. This fact, together with the approximate equality of 
the high-momentum tails of n{k) for all nuclei with A > 4, allows us to make the main assumption of 
this work namely, that the particle-state contributions to the momentum distribution are almost equal 
for all nuclei with A> 4. 

Let us decompose the proton momentum distribution (|^) in two terms: 



n{k) = nh{k) + np{k), 



(4) 



where the first term is the hole-state contribution 
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while the second one is the particle-state contribution 

^ oo 

np{k) = —Y^{2j + l)K Ra{k) 



(6) 



Using the assumed equality of the particle-state contributions np{k) for all nuclei, we obtain the following 
general relation of the correlated proton momentum distribution of a nucleus {A, Z) with that one of the 
'^He nucleus: 

Fa,. 
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n^'^{k) = N 
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where 



N = 



Fa.z 



nlj 



and Fa,z is the Fermi level for the nucleus (A,Z). 

Taking ^°Ca as an example the above expressions give: 
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As shown in , the hole-state natural orbitals are almost unaffected by the short-range correlations 
and, therefore, the functions Rnij{k) in Eq. (|^) can be replaced by the corresponding Hartree-Fock 
single-particle wave functions or by the shell- model single-particle wave functions Rnij{k). The hole- 
state occupation numbers Xnij are close to unity within the Jastrow correlation method ||2^ and we 
can set them equal to unity with good approximation. The properties of the hole-state natural orbitals 
and occupation numbers and the decomposition of the proton momentum distribution in the hole- and 
particle-state contributions (Eqs.(4)-(6)) lead to a similarity of our model to that one suggested for 
calculations of the spectral function in ||2^. In it the mean- field predictions for the spectral function 
are combined with its correlated part extracted from the nuclear matter calculations and recalculated 
for finite nuclei within the local density approximation. In our model, the correlated proton momentum 
distributions can be calculated for any nucleus by means of the occupied shell-model wave functions 
and the proton momentum distribution of the '^He nucleus which is taken from p^ ] and which contains 
short-range correlation effects. 



III. CALCULATIONS AND DISCUSSION 

In this work we calculate the proton momentmn distribution for nuclei ^^C, ^^O, '^^Ca, ^^Fe and ^^^Ph. 
Empirical estimations for n(k) are available for nuclei ^■^C and ^^Fe ||2^. 

In our calculations of proton momentum distributions we use two types of MFA single-particle wave 
functions: 1) single-particle wave functions obtained within the Hartree-Fock method by using Skyrme 
effective forces and 2) multiharmonic oscillator single-particle wave functions (with different values of the 
oscillator parameter for each state) which lead to a simultaneous description of ground-state radii and 
binding energies | |35| , p6[ |. In addition to p^ , in our calculations the multiharmonic oscillator s.p. wave 
functions arc orthonormalized. The values of all hole-state occupation probabilities Xnij in Eqs. (|^) and 
(H) are set equal to unity. The empirical data of n{k) for '^He are taken from j2^. They are given in 
Table I. As mentioned in the Introduction, the extraction of the data for n(k) is model-dependent. Due 
to this, we give in Fig.l the comparison of the data for n{k) in ^_ffe from [2^ with the calculations within 
the variational Monte Carlo method from [|l9|. As can be seen from Fig.l, the agreement is good and 
later we use in our calculations the data for n{k) in ^iJe from p6| . 
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The calculated proton momentum distributions for the nuclei examined are given in Figures 2-6, re- 
spectively. They are compared with the available data for ^^C and ^^Fe from |26) and the proton 
momentum distributions obtained in various theoretical methods, namely: for ^^C from p2| , for ^^O 
from |^,||,|9|-|3|], for "OCa from [|||3|] and for ^^^Pb from |3ll. For the ^^C nucleus the resuhs for 
the proton momentum distribution using the s.p. wave functions from the multiharmonic oscillator shell 
model but without including correlations are given in Fig. 7. Hence, the necessity of accounting for 
correlations becomes apparent. 

We have the following purposes within the practical method suggested in this work for realistic cal- 
culations of the nucleon momentum distribution in light, medium and heavy nuclei: 1) We like to show 
that the high-momentum tail of the momentum distribution for any nucleus can be approximated by 
that for ^He. We also check to what extent this approximation affects the central part of the momentum 
distribution. Since the low-momentum components of n{k) are determined mainly by the hole-state nat- 
ural orbitals contribution, the justification of the use of shell-model- or Hartree-Fock s.p. wave functions 
instead of hole-state natural orbitals can be checked; 2) We examine how well different s.p. shell-model 
wave functions can describe also the middle part of the momentum distribution which bridges the shell- 
model behaviour of the central part and the non-shell-model behaviour of the tail of the momentum 
distribution; 3) We like to apply this method in which correlation effects are extracted from *He to 
calculate n(k) as alternative one to the methods in which correlations are extracted from nuclear matter 
and in this way, if possible, to diminish the theoretical uncertainties on the momentum distribution for 
medium-heavy nuclei. 

One can see from Figs. 2 and 5 that the use of the single-particle wave functions from the multiharmonic 
oscillator shell model leads to better description of the experimental data for the central part of the 
momentum distribution than the use of the Hartree-Fock single-particle wave functions. In both cases 
the main deviations from the experimental data are for small momenta (fc < 0.5 fm^^). They are larger 
in the case when Hartree-Fock s.p. wave functions are used and this is a common feature of the results for 
all nuclei considered. This is due to the well-known fact that the Hartree-Fock method cannot give 
a realistic wave function for the Is state in the "^He nucleus. Namely this function {R^^'^ (k)) takes part 

1/2 

in the expression for n{k) (Eq.(7)) in all nuclei. Both types of s.p. wave functions, however, give similar 
results for the middle part as well as for the tail of the momentum distribution in all cases considered. 
The comparison of the results obtained by using of different mean-field single-particle wave functions 
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can be useful for the proper choice of the latter in the applications of the model to practical calculations 
of n(k) in cases when the knowledge of this quantity is necessary. 

Our numerical check shows that the results obtained by using values of the occupation probabilities 
Xnij coming either from the experiments or from nucleon-nucleon correlation methods (e.g. Jastrow one 
[p2[) are almost the same as those obtained with A„;j = l. As long as the correlated values of Xnij do not 
differ significantly from the value A„;j=l, the improvement is not sizeable. 

We emphasize that only hole-state occupation probabilities and wave functions enter the main relation- 
ships of the model (Eqs.(7) and (8)). In this way the suggested model can be easily applied to calculate 
momentum distributions in nuclei taking into account the nucleon-nucleon correlation effects. Concerning 
the particle-state contribution np{k) (Eq.(6)) to the proton momentum distribution (which is accounted 
for in the model by means of the term n ^'^(k) in Eq.(7)) we would like to mention that the decisive role 
for the existence of the high-momentum components in np{k) plays the form of the particle-state natu- 
ral orbitals Ra (which are strongly localized in coordinate space) but not the particle-state occupation 
numbers A„ {a > ap) 

As can be seen from Fig. 3 the results for the proton momentum distribution using the s.p. wave 
functions from the multiharmonic oscillator shell model are in agreement (at least for k > 0.8 fm^^) 
with those obtained in the variational Monte Carlo method ||2^ and in the calculations based on the 
local density approximation from p9[-|3l[| . The same can be seen from the comparison of our results for 
'^''Ca with those from given in Fig. 4. In our opinion, the calculations within the suggested model 
(with correlation effects from "^He) and the similarity of the results with those obtained in methods with 
correlations from nuclear matter can diminish in some sense the theoretical uncertainties on n{k). We 
would mention that some differences remain between our results and those from [3l| ] for ^^^Pb in the 
middle part of the momentum distribution (Fig. 6). 

It would be useful also to have estimations obtained extracting the correlation effects from nuclei 
ligther than ^He, such as the deuteron or ^He. As shown, however, in the variational correlation method 
(see Fig. 5 in pof), the high-momentum components of n{k) for deuteron and '^He, although with slopes 
similar with those of n{k) for '^He and nuclear matter, are quite different from them. At the same time, 
the high-momentum components of n{k) for ^He and nuclear matter are almost the same. This is the 
reason to use in our model correlation effects extracted from ^He and to relate them with the properties 
of the nuclear interior. 
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IV. CONCLUSIONS 



In the present paper a correlation model for calculating the proton momentum distribution in nuclei 
with ^ > 4 is proposed. The model combines the mean-field part of the momentum distribution with 
its correlated part taken from ^He on one and the same footing using the natural orbital representation. 
The estimation of the correlated part of n(k) is based on the well-known fact that the high-momentum 
components of the momentum distribution normalized to unity (at k > 2 fm^^) are nearly the same for 
all nuclei with A > 4. This fact, together with the use of the natural orbital representation, gives the 
possibility to obtain realistic momentum distributions in nuclei (including the regions of small momenta, 
k < 2 fm^^, and of the intermediate momenta, k ^ 2 fm^^) using only hole-state natural orbitals. 
The latter are replaced to a good approximation by shell- model single-particle wave functions. Thus 
the model gives a practical way for an easy calculation of the momentum distribution for any nucleus. 
The numerical results in this work confirm to a great extent the abilities of the suggested correlation 
model to give realistic estimations for the proton momentum distribution in ^^C and ^^Fe and to predict 
the behaviour of n{k) in ^®0, ''"Ca and '^^^Ph nuclei. They are in agreement with the results for the 
proton momentum distribution in ^^O and ^"Ca obtained within other theoretical methods in which the 
correlation effects are incorporated using nuclear matter results and with some empirical data for ^^C and 
^^Fe obtained using the y-scaling method. The knowledge of the realistic proton momentum distributions 
obtained in this work would allow us to describe in a similar way as it is done in quantities which 
are directly measurable in processes of particle scattering by nuclei. 
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Figure 1. Proton momentum distribution n{k) versus k of ^He. The solid triangles represent the data 
from The solid line is the result from The normalization is: J Ti{k)(fi'k = 1. 

Figure 2. Proton momentum distribution n(fc) versus k of ^^C. Calculations by using single-particle 
wave functions from the multiharmonic oscillator shell model p^ ] are presented by solid line and those 
by using Hartree-Fock single-particle wave functions by long-dashed line. The short-dashed line is n{k) 
calculated in the Jastrow correlation method |^^. The solid triangles represent the data from |^^. The 
normalization is as in Fig.l. 



Figure 3. Proton momentum distribution n(k) versus k of ^^O. The lines 1 and 2 are the results of 



the present work using Hartree-Fock and multiharmonic oscillator shell model |36| single-particle wave 
functions, respectively. The lines 3, 4, 5, 6 and 7 are the results from |30|, ||2^, |^ and pOt , 
respectively. The normalization is as in Fig.l. 



Figure 4. Proton momentum distribution n(k) versus k of '^'^Ca. The solid and long-dashed lines are as 
in Figure 2. The dotted line is n{k) calculated in ||2^. The short-dashed line is the result from |3^. The 
normalization is as in Fig.l. 

Figure 5. Proton momentum distribution n{k) versus k of ^^Fe. Calculations by using single-particle 
wave functions from the multiharmonic oscillator shell model js^] are presented by solid line and those 
by using Hartree-Fock single-particle wave functions by long-dashed line. The solid triangles represent 
the empirical data from |2^ . The normalization is as in Fig.l. 

Figure 6. Proton momentum distribution n{k) versus k of ^^^Pb. Calculations by using single-particle 
wave functions from the multiharmonic oscillator shell model js^] are presented by solid line and those 
by using Hartree-Fock single-particle wave functions by long-dashed line. The results from pi] ] are given 
by short-dashed line. The normalization is as in Fig.l. 

Figure 7. Proton momentum distribution n{k) versus k of ^^C . Calculations by using single-particle wave 
functions from the multiharmonic oscillator shell model |36| without including correlations are presented 
by solid line. The solid triangles represent the empirical data from p^. The normalization is as in Fig.l. 
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TABLE I. Experimental values of the proton momentum distribution in "^He |g6j. The normalization is: 



k n(k) k n(k) 

MeV/c /m^ MeV/c fm^ 

50 0.757916 350 0.001065 

100 0.331900 400 0.000749 

150 0.119684 450 0.000615 

200 0.039692 500 0.000548 

250 0.011019 550 0.000380 

300 0.002855 
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